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ABSTRACT 

Galaxy-scale outflows, which are thought to provide the link connecting the central 
black hole to its host galaxy, are now starting to be observed. However, the physical 
origin of the mechanism driving the observed outflows, whether due to energy-driving 
or radiation-driving, is still debated; and in some cases, it is not clear whether the 
central source is an active galactic nucleus (AGN) or a nuclear starburst. Here we 
study the role of radiation pressure on dust in driving galactic-scale AGN outflows, 
and analyse the dynamics of the outflowing shell as a function of the underlying 
physical parameters. We show that high-velocity outflows (>1000 km/s) with large 
momentum flux (> lOL/c) can be obtained, by taking into account the effects of 
radiation trapping. In particular, the high observed values of the momentum boosts 
can be reproduced, provided that the shell is initially optically thick to the reprocessed 
infrared radiation. Alternatively, the inferred measurements of the momentum flux 
may be significantly biased by AGN variability. In this context, the observations of 
powerful outflows on kiloparsec scales, with no or weak signs of ongoing nuclear activity 
at the present time, could be re-interpreted as relics of past AGN episodes. 
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1 INTRODUCTION 

The energy and momentum released by accretion onto the 
central black hole are thought to couple to the surrounding 
medium through some form of active galactic nucleus (AGN) 
feedback. As the accretion energy clearly exceeds the binding 
energy of the galaxy bulge, AGN feedback must have a great 
impact on its host, and may explain the observed black hole- 
host galaxy correlations (Silk & Rees 1998; Fabian 1999; 
King 2003; Murray et al. 2005; Fabian 2012, and references 
therein). Feedback from the central black hole is also invoked 
to account for a number of observed galaxy properties, in 
particular the quenching of star formation leading to the 
emergence of red and dead ellipticals (Springel et al. 2005; 
Di Matteo et al. 2005; Croton et al. 2006; Bower et al. 2006). 
In this context, AGN outflows are thought to provide the 
physical connection linking the small scales (^pc) of the 
central black hole to the large scales (>kpc) of the host 
galaxy where star formation takes place. 

There is increasing observational evidence of galactic- 
scale outflows reported in the literature (Fischer et al. 
2010; Feruglio et al. 2010; Rupke & Veilleux 2011; Sturm 
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et al. 2011; Veilleux et al. 2013; Cicone et al. 2014). These 
high-velocity, large-scale outflows are detected in different 
forms, including ionised, neutral, and molecular gas phases. 
Molecular outflows may be of particular interest since the 
molecular gas represents the bulk of the interstellar medium 
from which new stars eventually form. The measured high 
mass outflow rates, which often exceed the star formation 
rates, also imply short gas depletion timescales (e.g. Sturm 
et al. 2011; Cicone et al. 2014). This suggests that massive 
outflows can potentially remove gas from the entire galaxy, 
and thus critically affect the evolution of the host. Indeed, 
the detection of such powerful outflows has been posited as 
a proof of AGN feedback in action. The main question is 
what drives the observed outflows, and how the energy and 
momentum released by the central black hole actually cou¬ 
ple to the interstellar medium of the host galaxy. In prin¬ 
ciple, AGN feedback can operate via a number of physi¬ 
cal mechanisms, such as jets, winds, and radiation pressure 
(Fabian 2012, and references therein). One way of driving 
large-scale outflows is via high-velocity winds that gener¬ 
ate shockwaves propagating into the ambient medium (e.g. 
King et al. 2011; Zubovas V King 2012; Faucher-Giguere 
V Quataert 2012); another possibility is by direct radia¬ 
tion pressure on dust (e.g. Fabian 1999; Murray et al. 2005; 
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Ishibashi & Fabian 2012, see the latter for a discussion com¬ 
paring the two driving mechanisms). Focusing on radiative 
feedback, we note that the importance of radiation pressure 
on dust has already been highlighted in a number of past 
works. The direct effects of radiation pressure on dusty gas 
were first considered and introduced in the context of AGN 
feedback by Fabian (1999). Murray et al. (2005) discussed 
the role of radiation pressure on dust in driving galactic- 
scale winds through momentum deposition; while Thompson 
et al. (2005) examined the physics of starburst disks sup¬ 
ported by radiation pressure on dust grains. Numerical sim¬ 
ulations of radiative feedback, explicitly based on radiation 
pressure on dust, have also been carried out (Debuhr et al. 
2011; Roth et al. 2012; Krumholz & Thompson 2013; Davis 
et al. 2014). Most recently, Thompson et al. (2015) studied 
the dynamics of radiation pressure-driven shells and subse¬ 
quent applications to a wide variety of astrophysical sources. 

Here we wish to extend the work of Thompson et al. 
(2015) by further investigating the role of radiation pres¬ 
sure on dust in driving galactic-scale outflows, and discuss 
its implications in the broader context of AGN feedback. 
The paper is structured as follows. In Section 2, we recall 
some generalities on radiation pressure on dust as a poten¬ 
tial mechanism for driving large-scale outflows. In Section 
3, we study the case of the fixed-mass shell, and analyse the 
dependence of the shell dynamics on the different physical 
parameters. We then consider the sweeping-up of ambient 
material (Section 4) and the possibility of AGN variability 
(Section 5). We compare our model results with observa¬ 
tions of AGN-driven outflows in Section 6, and discuss the 
resulting implications in Section 7. 


2 RADIATION PRESSURE ON DUST 

The strength of the interaction between the AGN radia¬ 
tion and the surrounding medium depends on the underlying 
coupling process: electron scattering or dust absorption. It 
has been argued that transport by direct radiation pressure 
on electrons is problematic, as the galaxy becomes rapidly 
optically thin to the central radiation, and is indeed inade¬ 
quate to drive large-scale outflows. In contrast, the coupling 
can be much stronger in the case of radiation pressure acting 
on dust grains, as the dust absorption cross section (cr^) is 
much larger than the Thomson cross section (cr^), with a ra¬ 
tio typically of the order of ctc^/ctt ~ 1000 (e.g. Fabian 2012). 
Gorrespondingly, the effective Eddington limit for dust dif¬ 
fers from the standard Eddington limit by a factor cr^/crT. 

A major fraction of the AGN bolometric luminosity is 
emitted in the blue bump component, peaking in the ul¬ 
traviolet (UV) region. The UV photons are efficiently ab¬ 
sorbed by dust grains embedded in the gas, the typical UV 
opacity being of the order of Hiuv ~ lO^cm^g”^. As re¬ 
quired by energy conservation, the absorbed UV emission 
is then reprocessed and re-emitted as infrared (IR) radia¬ 
tion. Note that the typical opacity of dust in the infrared 
{i^iR ~ 5cm^g“^), although greater than the electron scat¬ 
tering opacity {rt — cftITU p ~ 0.4cm^g“^), is much smaller 
than at UV wavelengths (Novak et al. 2012). The re-emitted 
IR photons can then either directly escape the medium with¬ 
out further interactions (if the medium is optically thin), or 


undergo multiple scatterings before diffusing out of the re¬ 
gion (if the medium is optically thick). 

The momentum flux transferred from the radiation field 
to the gas is L/c in the single scattering limit. Thus each 
absorption and scattering of UV radiation couples L/c to 
the ambient medium. If the gas is optically thick to the 
re-radiated IR emission, additional momentum will be de¬ 
posited due to the multiple scatterings: TiRLjc where tir is 
the IR optical depth. 


2.1 Shell evolution 


We assume that radiation pressure drives a shell of dusty gas 
into the galaxy, and consider a thin shell geometry following 
Thompson et al. (2015). The general form of the equation 
of motion of the shell is given by: 


d 

dt 


^{Mshir)v) = ^(1 + TiR - e 


GM{r)Msh{r) 


( 1 ) 


where M{r) is the total mass distribution and Msh{r) is 
the shell mass. The total mass enclosed within a radius r is 
modelled as an isothermal sphere: 


= (2) 

where a is the velocity dispersion, and the additional term 
includes the self-gravity of the shell. The driving term is 
provided by the central luminosity L; the IR and UV optical 
depths are respectively given by 


riR(r) = 


RiRMshir) 

47rr^ 


Tuv{r) 


RuvMshjr) 

47rr^ 


(3) 

(4) 


where rir = 5 cm^g“^ fdg,MW and ruv = 10^cm^g“^ 
fdg,MW are the IR and UV opacities, with the dust-to-gas 
ratio normalized to the Milky Way value. 

As discussed in e.g. Thompson et al. (2015), there are 
three main regimes of optical depth to be identified: optically 
thick to both UV and IR, optically thick to UV but optically 
thin to IR, and optically thin to UV. If the medium is opti¬ 
cally thick to the reprocessed IR radiation, then most of the 
photons momentum can be transferred to the gas through 
radiative diffusion. Such conditions can be reached in the 
nuclear regions of dense starbursts and obscured AGN, typ¬ 
ically on scales of ~few hundred parsecs (Thompson et al. 
2005). In less extreme environments, the re-emitted IR pho¬ 
tons tend to leave the system without further interactions. 
In the single scattering regime (optically thick to UV but 
optically thin to IR), UV photons are absorbed, while the 
re-radiated IR photons freely escape. Such situations are 
met for a wide range of typical galaxies interstellar media. 
Eurther out, the ambient medium finally becomes optically 
thin to the central UV radiation. Due to their different dust 
opacities in the IR and UV regions, the IR transparency 
radius is always smaller than the UV transparency radius. 


3 FIXED-MASS SHELL 

We first consider the case of a geometrically thin shell of 
constant mass [Mshir) — Msh), starting from an initial ra¬ 
dius Ro with initial velocity vq , and moving outwards in the 
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galactic potential. In this simple case, analytical approxi¬ 
mations can be obtained (cf Thompson et ah 2015), which 
allow us to gain some physical insight into the problem. In 
order to study how the evolution of the shell depends on 
the different physical parameters, we consider the variation 
of one single quantity at a time (while keeping the other 
values fixed), and analyse its effects on the shell dynamics. 

We assume the following values as fiducial parameters 
of the model: a = 200 km/s, hur = 5 cm^/g, ruv = 10^ 
cm^/g, Ro = 200 pc, = 0 km/s, L = 5 x 10^®erg/s, 
Msh = 5 X lO^M©. For a shell of constant mass, the op¬ 
tical depth rapidly falls off with increasing radius, scal¬ 
ing as T oc 1/r^. The corresponding transparency radius 
for IR and UV radiation are respectively given by Rir = 

~ and Ruv = 


3 ^sh\. 7 ^P^ fiducial parameters. We note that the 
locations of the critical radii depend on the medium opacity 
and shell mass, scaling as Rir,uv oc . 


3.1 Radiative versus gravitational forces 


The evolution of the shell is governed by the balance between 
the inward force due to gravity and the outward force due 
to radiation pressure. The gravitational force is given by 


Fgrav ij ') 


GM{r)Msh 1 

--oc - 


( 5 ) 


and roughly scales as Fgrav oc 1/r for an isothermal mass 
distribution. The total radiative force is given by 


Fradij ') 


c 


( 6 ) 


The shell is initially optically thick to the re-radiated 
IR emission up to Rir^ beyond which it becomes optically 
thin. At small radii, the radiative force is thus dominated 
by a combination of the infrared and single scattering terms 

Frad{r) ^ ^{I+tir) (7) 

At larger radii, where the shell becomes optically thin 
to the UV radiation, the radiative force is dominated by the 
ultraviolet term 


Frad ('^) 


L 

~ '^UV — oc 
c 




( 8 ) 


scaling as 1/r^, and we recover the classical optically thin 
limit. 

In Figure 1, we plot the forces due to gravity (black 
curve) and radiation pressure (red curve) as a function of 
radius. We see that the radiative force exceeds the gravita¬ 
tional force for most radii, i.e. the ratio Frad jFgrav is greater 
than unity, implying that the shell is able to travel outwards. 
On large scales (r > Ruv), the radiative flux decreases more 
rapidly with radius (oc 1/r^) than gravity (oc 1/r), which 
eventually may dominate at very large radii. 



Figure 1. Gravitational force (black) and radiative force (red) 
as a function of radius. 



Figure 2. Critical luminosity as a function of radius. 


r 

which can be considered as a generalised form of the 
Eddington luminosity. The critical luminosity is plotted as 
a function of radius in Figure 2. We observe that LE{r) 
drops to a minimum around a few ~kpc, corresponding to 
the location of the maximal gap between the radiative and 
gravitational forces (i.e. maximum of the Frad j Fgrav ratio). 


3.1.1 Critical luminosity 

Equating the outward force due to radiation pressure to the 
inward force due to gravity, Frad{r) = Fgravir)^ leads to a 
critical luminosity: 


3.2 Asymptotic velocity 

Neglecting the exponential term in the equation of motion 
(Eq. 1), we obtain an analytical solution for the radial veloc¬ 
ity profile (cf Eq. 16 in Thompson et al. 2015). The resulting 
asymptotic velocity can be approximated as: 
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Figure 3. Velocity as a function of radius for variations in lu¬ 
minosity: L = 3 X 10"^® erg/s (blue), L = 5 x lO'^^erg/s (black), 
L = lO'^^erg/s (green). 



Figure 4. Velocity as a function of radius for variations in shell 
mass: Mgh = 10 ® Mq (blue), Mgh = 3 x 10 ® Mq (black), Mgh = 
10 ®M 0 (green). 


HiUV ^ 1/4 


c ^AirMsh 


) ^ OC Kjjy 




( 10 ) 


Numerically integrating the equation of motion, we ob¬ 
tain the exact radial velocity profile of the outflowing shell. 
We note that the isothermal potential approximation only 
holds in the inner regions (on ~10 kpc-scales), but the re¬ 
sults remain very similar assuming a NFW profile. Thus in 
the following, we show for simplicity the results obtained in 
the case of the isothermal potential. Figures 3 and 4 show 
the velocity of the shell as a function of radius for different 
values of the central luminosity and shell mass, respectively. 
We see that the shell reaches higher velocities for higher 
luminosities, as expected. For a given luminosity, higher ve¬ 
locities are reached by lower mass shells. A variation in the 
central luminosity has a stronger effect on the velocity pro¬ 
file than a change in the shell mass. As can be seen from 
the analytical limit (Eq. 10), we note that the dependence 
of the asymptotic velocity on central luminosity is stronger 
[v OC L^/^) than on shell mass [v oc For a given cen¬ 

tral luminosity and shell mass, the asymptotic velocity also 
depends on the UV opacity of the medium. Larger values of 
nuv lead to higher velocities, with considerable differences 
seen at large radii; an enhanced IR opacity also gives a small 
contribution to the increase in shell velocity. 

Considering variations in the shell initial conditions, we 
observe that the initial velocity has no significant effect on 
the shell dynamics, while a smaller initial radius leads to a 
slightly higher asymptotic velocity. Similarly, the shell at¬ 
tains higher velocities in systems with lower velocity dis¬ 
persions, as might be expected in shallower gravitational 
potentials. 


3.3 Momentum ratio 

The momentum ratio of the shell is dehned as 


c = 


Msh'^ 


( 11 ) 


At large radii (r > Ruv), the asymptotic momentum 
ratio is given by 

C(r > Ruv) ~ 2^ ~ (12) 

r V TT r 

and scales as 1/r. 

At small radii, close to the IR transparency radius (r ~ 
Rir), the momentum ratio can be approximated as 

(.3, 

The maximal value of the momentum hux is attained 
at small radii (r < Rir), where the shell is optically thick 
to the infrared radiation. If the initial radius lies close to 
the IR transparency radius (Rq ~ Rir), then the shell be¬ 
comes rapidly optically thin to the re-radiated IR emission 
and C{Rir) ~ 2. We note that the momentum ratio can 
be signihcantly higher only if the optical depth to the re¬ 
radiated IR photons is much larger than unity at the launch 
radius. As we will see later, the IR optical depth due to radi¬ 
ation trapping is indeed the main parameter governing the 
momentum ratio of the outhowing shell. 

Figure 5 shows the momentum ratio of the shell as a 
function of radius for different values of the central lumi¬ 
nosity. We see that the momentum ratio always decreases 
as 1/r at large radii, while the peak values are roughly of 
the same order {(max ~ 1). A qualitatively similar plot is 
obtained in the case of varying shell mass. Thus, in contrast 
to the previously discussed velocity prohles, the momentum 
ratio is not very sensitive to variations in central luminosity 
and/or shell mass. 

On the other hand, variations in the initial IR optical 
depth can have a signihcant impact on the momentum hux 
of the outhowing shell, with a larger initial IR optical depth 
leading to signihcantly higher ^ values. In physical terms, 
the increase in the initial IR optical depth can be due either 
to an enhanced IR opacity or a smaller initial radius. An 
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radius (kpc) 


Figure 5. Momentum ratio as a function of radius for variations 
in luminosity: L = 3x lO'^^erg/s (blue), L = 10^®erg/s (black), 
L = lO'^^erg/s (green). 



Figure 6. Momentum ratio as a function of radius for variations 
in initial radius: Rq = 200 pc (black), Rq = 50 pc (blue), Rq = 
20 pc (green). 


enhanced opacity can be obtained e.g. by assuming a higher 
dust-to-gas ratio, as might be the case in the central regions 
of dense ULIRGs. In Figure 6 we plot the radial profile of 
the momentum ratio for 3 different values of Rq . We see that 
larger momentum ratios are obtained for smaller initial radii. 
We also note that high ( values are only observed on small 
scales (<kpc); while at larger radii (>kpc), the momentum 
flux profiles rapidly fall off as oc 1/r. From Eq. 13, we note 
that C(^/-r) scales inversely with the initial radius Rq. A 
lower limit on the initial radius is set by the dust sublimation 
radius given by 


■ / IGnasBT^^, 

where Tgub is the sublimation temperature of the dust. For 
typical values (Ts^ife=1500K), the sublimation radius is of 
the order of Rsub ~ Ipc. 

We also verify that variations in the other parameters 
do not affect the momentum ratio of the shell in a significant 
way. In particular, changes in the UV opacity do not have 
an appreciable effect on the value of the momentum flux. 

Summarising, a favourable combination of the different 
physical parameters can lead to high velocity outflows with 
large momentum fluxes. For a given luminosity and shell 
mass, the asymptotic velocity of the shell is mainly deter¬ 
mined by the UV opacity; whereas the momentum flux (and 
in particular its peak value) is most sensitive to variations 
in the IR optical depth. On the other hand, we note that 
changes in the UV opacity have no significant impact on 
the momentum ratio; conversely, variations in the IR op¬ 
tical depth do not much affect the asymptotic velocity of 
the shell. These results are in qualitative agreement with 
the analytical limits given in Eqs. 12 and 13. We can thus 
disentangle the respective contributions of the UV and IR 
terms to the overall dynamics of the outflowing shell. 



Figure 7. Mass outflow rate as a function of radius for different 
shell models: flxed mass shell with = 5 x 10® Mq (black), 

shell sweeping up mass from a constant density medium with 
no = 0.1cm“® (blue), shell sweeping up mass from an isothermal 
density distribution with fg = 0.05 (green). 


3.4 Outflow rate 


The mass outflow rate in the shell limit can be defined as: 



MshV 

r 


(15) 


In Fig. 7, we plot the mass outflow rates as a function 
of radius for different shell models (cf Sect. 4). In the case 
of a fixed mass shell, the outflow rate can reach values of 
several hundred solar masses per year, followed by a decline 
scaling as M oc 1/r (black curve). 
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4 EXPANDING SHELLS 

In more realistic situations, the shell is likely to sweep up 
mass as it expands into the interstellar medium of the host 
galaxy. The previous analysis of the dynamics of hxed-mass 
shells may be a good approximation in cases where the sub¬ 
sequently added swept-up mass (Msw{r)) is negligible com¬ 
pared to the initial mass of the shell (Msh.o)- But when the 
former becomes comparable to or exceed the initial mass, 
one has to explicitly take into account the contribution of 
the swept-up material. More specifically, the mass of the 
swept-up material on scales of the UV transparency radius 
{Msh{Ruv)) should be compared with the initial mass of 
the shell (as discussed in Thompson et al. (2015)). 

The density distribution of the ambient medium can be 
parametrized as a power law of radius with slope a: 

n{r) = no (^^) (16) 



where no is the density of the external medium. The corre¬ 
sponding swept-up mass is given by: 


Tfstu (^ ) — ^iTlTTlp 



dr = AirmpnoRQ 


3 — a 


(17) 


In the following, we focus on two particular cases of the 
ambient density profiles: a = 0 and a = 2, corresponding to 
constant density and isothermal distributions, respectively. 
The dust-to-gas ratio of the swept-up matter is parametrized 
as ^ = fdg/fdg,MW, normalized to the Milky Way value. 


4.1 Constant density distribution 

If the ambient medium is characterized by a constant 
density distribution {a = 0), the corresponding swept- 
up mass is given by Msw,cst{r) = ^Tvmpnor^ = 2 x 
10^no, 2 rlpcMQ. The critical density for which the swept-up 
mass on scales of Ruv is equal to the initial shell mass is: 
Uc = . As the swept-up mass increases 

very rapidly with radius (due to the steep radial dependence 
of the form Msw,cst{r) oc r^), it quickly becomes dominant 
at large radii, even for external densities lower than the crit¬ 
ical value (cf Thompson et al. (2015)). As in the case of the 
fixed-mass shell, variations in the physical parameters lead 
to qualitatively similar effects on the shell dynamics. For 
instance, outflowing shells reach higher velocities with in¬ 
creasing luminosity and decreasing mass. 

Here we focus on variations in the value of the external 
density (no). From the plot of the radial velocity profile (Fig. 
8), we see that the shell reaches a certain maximal speed 
and then slows down, due to the sweeping up of ambient 
material. At small radii, the swept-up mass is still negligible 
compared to the initial mass of the shell, and the v(r) curves 
are roughly identical. In fact, the Msh{r)/Msh,o mass ratio 
starts to exceed unity at relatively large radii, around Vc = 
8 . 7 ^P^’ scales comparable to Ruv- 

The peak velocity attained and the subsequent decel¬ 
eration pattern depend on the density of the surrounding 
medium. As the external density increases, the shell sweeps 
up more mass and attains a lower peak velocity followed by 
a rapid deceleration. Since the velocity rapidly falls off at 
large radii, it is unlikely that the shell is able to escape the 
galactic halo, and plausibly falls back. On the other hand. 


Figure 8. Velocity as a function of radius for various values of 
spatially constant density: no = 0.5cm“^ (black), no = 0.1cm“^ 
(blue), no = 0.05cm“^ (green). 



Figure 9. Momentum ratio as a function of radius for various 
values of spatially constant density: no = 0.5cni“^ (black), no = 
0.1cm“^ (blue), no = 0.05cm“^ (green). 


we observe that variations in the values of the external den¬ 
sity do not have a dramatic effect on the momentum ratio 
of the shell (Fig. 9). 

In Fig. 7, we also plot the outflow rate for a shell sweep¬ 
ing up material from the surrounding medium of constant 
density (blue curve). As expected, the outflow rate scales as 
M (X r^v ai large radii, where the swept-up mass dominates 
the shell dynamics. 


4.2 Isothermal density distribution 

If the ambient density profile follows an isothermal distri¬ 
bution (a = 2), the swept-mass is given by Msw,iso{r) ~ 


© 2012 RAS, MNRAS 000 , 1-11 












AGN radiation pressure-driven outflows 7 


r = 2 X 10^ fg,o.oicr2oorkpcMQ, where fg is the gas 
fraction. The critical value of the gas fraction for which 
the swept-up mass at Ruv equals the initial mass is: fc = 
y. Again, varying the different physical 
parameters lead to qualitatively similar outcomes to the case 
of fixed-mass shells. 

We concentrate here on variations in the gas fraction, 
which is the most important parameter in determining the 
shell dynamics. As expected, the shell velocity is lower for 
ambient medium having higher gas fractions (Fig. 10). For 
fg > 0.05, the shell can never escape the galaxy as its veloc¬ 
ity is always lower than the escape velocity (v < 2a). For a 
smaller gas fraction {fg = 0.01), the shell reaches a higher 
peak velocity temporarily exceeding the escape speed, but 
later falls back. This suggests that the outflow is bound on 
large scales, although the actual escape condition may de¬ 
pend on the large-scale gravitational potential. In contrast 
to the velocity profiles, we note that the momentum ratio 
is not significantly affected by variations in the gas fraction 
(Fig. 11). Compared to the constant density case, the swept- 
up mass in the isothermal density case rises more gently with 
radius, following a linear relation of the form Msw,iso{r) oc r. 
The amount of swept-up material in the two cases are com¬ 
parable over a limited range of radii around Ruv, for the 
spanned external densities and gas fractions. The swept-up 
masses for the constant and isothermal density cases quickly 
diverge at larger radii, due to their distinct radial depen¬ 
dence. In Fig. 7, we also plot the outflow rate for a shell 
sweeping up matter from the surrounding environment with 
an isothermal density distribution (green curve). We observe 
that the mass outflow rate is proportional to the velocity 
(M oc t;), and is of the order of a few hundred solar masses 
per year. 

The decrease of the gas mass fraction over cosmic time 
has been reported in a number of works (Daddi et al. 2010; 
Tacconi et al. 2010). High-redshift galaxies tend to have gas 
fractions much higher than the values observed in present- 
day galaxies. The depletion of gas could be partly due to 
past outflow events, which may lead to lower gas fractions 
at later times, which in turn may facilitate subsequent out¬ 
flow episodes. In fact, cool gas may be removed from the 
core of the galaxy, leaving the central parts temporarily 
depleted. But the outflowing gas is likely to be bound on 
larger scales and remain trapped in the halo and surround¬ 
ing circum-galactic medium. Indeed, recent observations in¬ 
dicate the presence of significant amounts of cool gas sur¬ 
rounding quasar host haloes (Prochaska et al. 2013). Several 
gas re-accretion events may be expected over cosmic time, 
with a complex combination of escape and fail-back episodes 
determining the subsequent evolution. 


5 AGN VARIABILITY: TIME DEPENDENCE 

The observed outflow characteristics may also be affected by 
the temporal variability of the central AGN. Indeed, AGNs 
are known to vary by several orders of magnitude in lumi¬ 
nosity over a wide range of timescales. Radiation pressure 
on dust can potentially disrupt the ambient material, such 
that the fuel supply for the central quasar may be blown 
away following an outburst event. The temporal variability 
of the luminosity output can have important consequences 



Figure 10. Velocity as a function of radius for various values of 
the gas fraction: fg = 0.1 (black), fg = 0.05 (blue), fg = 0.01 
(green). 



Figure 11. Momentum ratio as a function of radius for various 
values of the gas fraction: fg = 0.1 (black), fg = 0.05 (blue), 
fg = 0.01 (green). 


on the inferred value of the momentum flux that should be 
taken into account. 

Observationally, the outflow is detected and its velocity 
measured at a certain distance from the centre, while the 
central luminosity is derived from the current observed flux. 
These two quantities are then used to compute the momen¬ 
tum ratio defined in Eq. 11. In order for the inferred values of 
C to be accurate, one requires the central luminosity to stay 
constant for a few million years, which is quite unlikely due 
to the variable nature of the AGN. Indeed, a certain amount 
of time is needed for the shell to reach the current location, 
and over this same time span the central luminosity may 
not have remained the same as when the shell was initially 
launched. If e.g. the central AGN luminosity has consider- 
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ably decreased over this time lapse, then the measured value 
of the momentum flux may be largely over-estimated. 

For simplicity, we hrst consider the case in which the 
central AGN undergoes a sudden drop in luminosity at a 
given time (Fig. 12). Despite the luminosity decrease, the 
outflow still carries on, but an abrupt jump occurs in the 
the corresponding value of the momentum flux. There are 
two main parameters governing the momentum ratio prohle: 
the amplitude of the luminosity drop and the time when this 
break occurs. If the subsequent luminosity (L 2 ) is lower, then 
the peak value of the momentum ratio is correspondingly 
higher; while for a shorter tb, the break occurs at a smaller 
radius. Thus, varying the two parameters, one can obtain 
different peak locations and distinct peak values. 

We next consider the case of an exponential decay in 

_ ^ 

luminosity: L(t) = Loe , where Lq is the initial luminosity 
and tc is the characteristic timescale (Fig. 13). We see that 
relatively large values of the momentum flux (C ^ 1) can 
be reached on small scales, provided that the characteristic 
timescales are quite short, of the order of a few million years 
{tc ~ lO^yr). Considering that AGN accretion is limited by 
self-gravity (cf Discussion), the characteristic timescale is 
mainly determined by the mass enclosed within the self¬ 
gravity radius, and tc may be of the order of ~ lO^yr at the 
Eddington limit (King & Pringle 2007). 

The time dependence of the AGN luminosity output 
may be relevant in interpreting the observational results, in 
particular for sources accreting at sub-Eddington rates at 
the present time, which may form a significant fraction of 
the galaxy samples discussed in the literature (Veilleux et al. 
2013; Cicone et al. 2014). In fact, there are several observa¬ 
tional indications which suggest variations by a few orders of 
magnitude in luminosity over relatively short timescales of 
~ (10^ — 10®)yr (Hickox et al. 2014, and references therein). 
A nice example is given by the so-called Hanny’s Voorwerp 
system, in which the central AGN luminosity dropped by 
a factor of > 10^ over the last ~ lO^yr, leaving the high- 
ionization cloud as a unique remnant of past activity (Keel 
et al. 2012). 


6 COMPARISON WITH OBSERVED AGN 

OUTFLOWS 

A number of detections of massive outflows extending on 
galactic scales have been reported over the past few years. 
Observations of molecular outflows, tracing the bulk of the 
interstellar medium from which new stars form, are of par¬ 
ticular relevance in determining whether the central AGN 
plays a role in quenching star formation. 

One of the best studied cases is Mrk 231, a nearby 
quasar in late stage merger, which is one of the most lu¬ 
minous ULIRGs in the local Universe (Eischer et al. 2010; 
Eeruglio et al. 2010; Rupke & Veilleux 2011). Observations 
of OH features and GO lines indicate the presence of a mas¬ 
sive molecular outflow, spatially resolved on ~kpc scale and 
with an outflow rate largely exceeding the star formation 
rate (Eischer et al. 2010; Eeruglio et al. 2010). A powerful 
neutral outflow, with velocities reaching ~ llOOkm/s and 
extending out to radii of 2-3 kpc, is also detected in this 
object (Rupke V Veilleux 2011). 

Similarly, Herschel-PAGS observations of local ULIRGs 



Figure 12. Momentum ratio as a function of radius in the case 
of a sudden drop in luminosity: Li = 5 X 10^®erg/s, L 2 = 1 x 
1046erg/s, tb = 5 X lO^yr (black); Li = 5 x 10^®erg/s, L 2 = 
5 X lO'^^erg/s, tb = 2 x lO^yr (magenta). 



Figure 13. Momentum ratio as a function of radius in the case of 
an exponential decay in luminosity: Lq = lO'^^erg/s, G = lO^yr 
(black), tc = ^ X lO^yr (blue), G = lO^yr (green). 


reveal the presence of massive molecular outflows with ve¬ 
locities reaching v ~ lOOOkm/s and outflow rates up to 
M ~ lOOOMo/yr (Sturm et al. 2011). The detection of 
molecular outflows has further been conhrmed in a larger 
representative sample of local galaxy merger systems, in¬ 
cluding ULIRG and QSO (Veilleux et al. 2013). Recent re¬ 
sults based on GO observations indicate that such powerful 
molecular outflows commonly have high momentum fluxes, 
of the order of ~ 20L/c (Gicone et al. 2014). However, one 
should recall that the quoted values are based on uniform 
spherical geometry, while in the shell limit the corresponding 
values should be divided by a factor of 3 (Gicone et al. 2014; 
Maiolino et al. 2012). Moreover, one should keep in mind 
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that the observational samples are biased toward sources 
with previously known outflow detections. 

In our picture, large-scale, high-velocity outflows may 
be obtained by considering the effects of radiation pressure 
on dust. We have seen that radiation pressure-driven shells 
can reach asymptotic velocities of ~1000 km/s on kpc-scales, 
mainly depending on the central luminosity and shell mass. 
In particular, a higher asymptotic velocity is reached for 
low-mass shells in high luminosity systems, with a stronger 
dependence on the central luminosity. Sturm et al. (2011) 
hnd a tentative evidence for a positive correlation between 
terminal velocity and AGN luminosity. A similar trend of 
increasing outflow velocity with increasing AGN luminosity 
is also reported in the larger sample of Veilleux et al. (2013), 
a trend in qualitative agreement with our results. 

On the other hand, in order to account for the high 
observed values of the momentum boosts, a number of 
favourable conditions need to be met. Following our param¬ 
eter space study, the main factor governing the small-scale 
dynamics is the IR optical depth. We have seen that a high 
initial IR optical depth, either due to an enhanced IR opacity 
or a smaller initial radius, can lead to large momentum ra¬ 
tios. In fact, the highest C values are attained at small radii, 
where the medium is optically thick to the reprocessed IR 
emission, i.e. in the obscured phase. For a given IR opacity, 
the momentum flux is most significantly affected by varia¬ 
tions in the initial radius: a decrease in Rq can lead to a 
significant increase in on ^kpc scales (cf Fig. 6). In this 
way, one can obtain large values of the momentum flux by 
adjusting the initial radius, without requiring extreme val¬ 
ues of the IR opacity. Based on the observed source proper¬ 
ties and inferred measurement of the momentum flux, one 
can in principle derive the initial radius. As a specific ex¬ 
ample, we can consider the case of Mrk 231. Adopting the 
reported values for the velocity dispersion (a = 120km/s), 
central luminosity (L ~ 10^®erg/s), and outflowing mass 
{Msh ~ 3x 10^Mq) (Veilleux et al. 2013; Gicone et al. 2014), 
we obtain values of C ^ 10 at the radius of observation, con¬ 
sistent with the value of = 26/3 reported in Gicone et al. 
(2014), for an initial radius of Rq < lOpc. Thus the mo¬ 
mentum flux may provide a constraint on the initial radius, 
requiring a quite small value of Rq in this particular system 
(for a hducial IR opacity). 

In principle, such method can be applied to other 
sources and used to infer the underlying physical quantities 
from the observed outflow parameters. However, the imple¬ 
mentation may not be straightforward due to the uncertain¬ 
ties and degeneracies involved, and especially in cases where 
the AGN luminosity has considerably dropped since the time 
when the outflow was initially launched. This may be par¬ 
ticularly true in systems currently accreting at very sub- 
Eddington rates, which may comprise a significant fraction 
of the observational samples (Veilleux et al. 2013; Gicone 
et al. 2014). Veilleux et al. (2013) estimate that sources in 
their sample accrete at sub-Eddington rates, in the range 
LjLE ~ 0.04 — 0.4. Similarly, the majority of the objects 
discussed in Gicone et al. (2014) have rather low Eddington 
ratios, with some sources accreting at very sub-Eddington 
rates. Eor instance, Mrk 273 with an estimated black hole 
mass of Mbh ^ 10^M© (U et al. 2013) should now be accret¬ 
ing below 1% of its Eddington limit. As discussed in Section 
5, variations in the central AGN luminosity can seriously af¬ 


fect the inferred measurement of the momentum boost. Since 
large drops in luminosity are commonly observed (Hickox 
et al. 2014, and references therein), correspondingly large 
values of C could be easily obtained. Observations report 
typical momentum rates of ~ lOL/c for molecular outflows 
extending on ~kpc scales (Gicone et al. 2014). In order to 
obtain high values of C at small radii, one requires short 
characteristic timescales (cf Eig. 12). 


7 DISCUSSION 

Galaxy-scale outflows, extending on ~kpc scales, are now 
commonly detected in a number of active galaxies (e.g. 
Sturm et al. 2011; Veilleux et al. 2013; Gicone et al. 2014). 
The observed molecular outflows are typically characterised 
by high velocities (v ~ lOOOkm/s), large mass outflow rates 
(M ~ lOOOMo/yr), and high momentum fluxes L/c). 
The inferred energetics of such powerful outflows suggest 
that they are most likely driven by the central AGN. But 
the problem remains as how to efficiently couple the energy 
and momentum released by the accreting black hole to the 
surrounding medium. 

As mentioned in the Introduction, large-scale outflows 
can be driven by a number of physical mechanisms. Eor 
instance, two distinct regimes of outflow driving are dis¬ 
cerned in the shocked wind scenario (e.g. King et al. 2011): 
momentum-driving at small radii (r < Ikpc), where the 
Gompton cooling time is shorter than the wind flow time, 
and energy-driving at larger radii where the full energy of 
the wind is coupled to the ambient medium. The energetics 
of the observed outflows, and in particular the high mea¬ 
sured values of the momentum ratios (~ 20L/c), seem to 
favour the energy driving mechanism (Zubovas V King 2012; 
Eaucher-Giguere V Quataert 2012) and apparently rule out 
driving by direct radiation pressure. 

However, radiation pressure on dust may be another vi¬ 
able mechanism for driving fast outflows on galactic scales. 
We note that radiation pressure acts on the bulk of the 
mass and do not expect developments of strong shocks in 
the medium. Radiative cooling of dense gas should also be 
efficient, allowing the dusty gas to cool down and remain 
cold. We have seen that shells driven by radiation pressure 
on dust can actually attain high velocities and large values of 
the momentum flux under certain physical conditions. The 
main parameter determining the small-scale physics and the 
peak of the momentum ratio is the initial IR optical depth. 
Indeed large values of the momentum flux consistent with 
the observed ones can be obtained, provided that the shell is 
initially optically thick to its own reprocessed IR radiation. 
In fact, high values of tjr are not unexpected in the cen¬ 
tral regions of dense ULIRG and obscured AGN (Thompson 
et al. 2005). 

In general, the optical depth depends on the density 
and opacity of the ambient medium. In realistic situations, 
the density distribution is likely to be inhomogeneous and 
clumping of the interstellar medium may lead to a lower¬ 
ing of the global dust covering fraction. As a consequence, 
radiative momentum transfer may be less efficient, and the 
optical depth might be reduced along certain lines of sight. 
Numerical simulations of radiative feedback indicate that 
the radiation force on dusty gas may be reduced by a factor 
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of ~ 2 in cases of severe clumping compared to the case of 
a smooth gas distribution (Roth et ah 2012). Although the 
initial launch radius may in principle be derived from the 
outflow parameters (Sect. 6), a clumpy distribution also im¬ 
plies that the initial radius is not strictly defined, and the 
overall picture may be much more subtle. A lower limit is 
only set by the dust sublimation radius, below which elec¬ 
tron scattering dominates the local opacity. 

The coupling between radiation and matter depends on 
the opacity of the ambient medium, which is directly pro¬ 
portional to the dust-to-gas ratio. An increase in the dust- 
to-gas fraction leads to an enhanced opacity, which in turn 
leads to a higher optical depth. Recent observations indicate 
that large amounts of dust can be produced in core-collapse 
supernovae and spread into the surrounding environment, 
implying high dust-to-gas ratios (Wesson et al. 2015; Owen 
Sz Barlow 2015). Moreover, their fits to the spectral energy 
distribution favour large dust grains, which are also more 
likely to survive. 

Hydrodynamic simulations of galaxy mergers with feed¬ 
back based on radiation pressure on dust require large val¬ 
ues of the IR optical depth (tir ~ 25) in order to repro¬ 
duce the observed M — a relation (Debuhr et al. 2011). This 
suggests that much of the black hole growth occurs in the 
regime where matter is optically thick to the re-radiated 
IR emission, i.e. in the obscured phase (cf Fabian 1999). 
Radiative transfer calculations, including multi-dimensional 
effects, indicate that winds driven by radiation pressure on 
dust can effectively attain momentum rates of several (~ 5) 
L/c due to multiple scatterings (Roth et al. 2012). On the 
other hand, Novak et al. (2012) note that the low opacity 
of dust to the re-radiated IR emission and dust destruction 
processes (e.g. sputtering) do not allow values much exceed¬ 
ing ~ L/c. It has also been argued that Rayleigh-Taylor 
instabilities developing in the outflowing medium tend to 
reduce the matter-radiation coupling and thus limit the mo¬ 
mentum transfer, such that the outflow cannot effectively 
exceed the single scattering limit (Krumholz Sz Thompson 
2013). By contrast, the most recent results show that, de¬ 
spite the development of instabilities in the flow, there is a 
continues acceleration of the gas, which potentially allows 
the propagation of the outflow on much larger scales (Davis 
et al. 2014). Therefore, although several issues still remain 
to be resolved, radiation pressure on dust should be consid¬ 
ered as a viable mechanism for driving large-scale outflows 
in AGN host galaxies. In fact, we have previously argued 
that the basic properties of galaxies, such as characteristic 
radius and mass, could be essentially set by the action of 
radiation pressure on dust (Ishibashi & Fabian 2014). 

Another concern is the variability of the central AGN. 
As we have seen in Section 5, luminosity variations can have 
considerable effects on the inferred value of the momentum 
ratio: if the current luminosity is much lower than the lu¬ 
minosity at the time when the shell was initially launched, 
then the momentum flux can be largely over-estimated. One 
may argue that AGN variability is an intrinsically random 
process, whereby luminosity variations occur in both direc¬ 
tions with essentially equal probability. But there could also 
be some preferential trend, with one way being more likely 
than the other. If the outflow sweeps up and carries away 
large amounts of accreting matter, then potential fuel is re¬ 
moved from the nuclear region. As a result, the accretion 


rate onto the central black hole may decrease, which in turn 
may lead to a lower luminosity output. Gertainly, the ac¬ 
cretion process is a complex phenomenon, and one needs 
to consider more in detail the connection with the outflow 
component. 

The AGN variability is also related to another inter¬ 
esting question regarding the outflow driving mechanism: 
AGN- or starburst-driven. When powerful outflows are ob¬ 
served on ~kpc-scales, but no corresponding sign of AGN 
activity is detected at the centre, one may be led to conclude 
that these outflows are not AGN-driven. Indeed observations 
of high-velocity outflows from (post)-starburst galaxies have 
been explained in terms of compact starbursts, without the 
need to invoke AGN feedback (Diamond-Stanic et al. 2012; 
Sell et al. 2014). However, the lack of ongoing AGN activity 
at the present epoch does not preclude that the large-scale 
outflows we observe today were not driven by a past pow¬ 
erful AGN episode that has since faded. In fact, there are 
now indications of the intermittent and episodic nature of 
AGN activity, with luminosities dropping by several orders 
of magnitude on relatively short timescales (Hickox et al. 
2014, and references therein). Recent ALMA observations of 
an obscured quasar, presenting both compact molecular and 
extended ionized outflows, have been interpreted in terms of 
episodic AGN activity whereby the two outflow components 
may be driven by two distinct burst events (Sun et al. 2014). 
An interesting case is the recent claim of an extended, high- 
velocity molecular outflow from a compact starburst (Geach 
et al. 2014). This has been attributed to stellar feedback, as 
the outflow momentum flux is found to be comparable to 
the radiative momentum output L/c from stars. However, if 
the bulk of the outflow only subtends a small solid angle at 
the starburst (as seems to be the case from the galaxy im¬ 
age, see their Fig. 3), then the comparison may break down, 
unless the starburst emission is particularly anisotropic. 

Given that the maximal luminosity of an AGN is likely 
within a factor of a few of the Eddington luminosity, a plau¬ 
sible limit on the momentum flux may be obtained by con¬ 
sidering the Eddington limit of the AGN rather than the es¬ 
timate of the current luminosity. The luminous AGN phase 
may last until the accretion disc empties. The actual ex¬ 
tent of the disc is likely limited by self-gravity, with gas 
beyond the self-gravity radius {Rsg) being consumed in star 
formation rather than feeding the central black hole (King 
Sz Pringle 2007). Such accretion episodes are shown to fol¬ 
low a characteristic time-evolution, with the disc evolution 
timescale {rgg) only weakly depending on the black hole 
mass (cf Eq. 14 in King & Pringle (2007)). 

In general, understanding the connection between the 
central AGN activity and star formation on galactic scales is 
a complex problem, and it is often not trivial to disentangle 
the underlying mechanism driving the observed outflows. In 
a number of situations, we may re-interpret the fast outflows 
attributed to central starbursts as being driven by past AGN 
activity. If this is indeed the case, then AGN feedback may 
play an even more important role in launching galactic-scale 
outflows than previously thought. 
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